Early-stage idiopathic Parkinson's disease is associated with reduced circular RNA expression.

Neurodegeneration in Parkinson's disease (PD) precedes diagnosis by years. Early neurodegeneration may be reflected in RNA levels and measurable as a biomarker. Here, we present the largest quantification of whole blood linear and circular RNAs (circRNA) in early-stage idiopathic PD, using RNA sequencing data from two cohorts (PPMI = 259 PD, 161 Controls; ICICLE-PD = 48 PD, 48 Controls). We identified a replicable increase in TMEM252 and LMNB1 gene expression in PD. We identified novel differences in the expression of circRNAs from ESYT2, BMS1P1 and CCDC9, and replicated trends of previously reported circRNAs. Overall, using circRNA as a diagnostic biomarker in PD did not show any clear improvement over linear RNA, minimising its potential clinical utility. More interestingly, we observed a general reduction in circRNA expression in both PD cohorts, accompanied by an increase in RNASEL expression. This imbalance implicates the activation of an innate antiviral immune response and suggests a previously unknown aspect of circRNA regulation in PD.

Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

Mismatch repair deficiency testing in Lynch syndrome-associated urothelial tumors.

Introduction: Lynch syndrome-associated cancer develops due to germline pathogenic variants in one of the mismatch repair (MMR) genes, MLH1, MSH2, MSH6 or PMS2. Somatic second hits in tumors cause MMR deficiency, testing for which is used to screen for Lynch syndrome in colorectal cancer and to guide selection for immunotherapy. Both MMR protein immunohistochemistry and microsatellite instability (MSI) analysis can be used. However, concordance between methods can vary for different tumor types. Therefore, we aimed to compare methods of MMR deficiency testing in Lynch syndrome-associated urothelial cancers. Methods: Ninety-seven urothelial (61 upper tract and 28 bladder) tumors diagnosed from 1980 to 2017 in carriers of Lynch syndrome-associated pathogenic MMR variants and their first-degree relatives (FDR) were analyzed by MMR protein immunohistochemistry, the MSI Analysis System v1.2 (Promega), and an amplicon sequencing-based MSI assay. Two sets of MSI markers were used in sequencing-based MSI analysis: a panel of 24 and 54 markers developed for colorectal cancer and blood MSI analysis, respectively. Results: Among the 97 urothelial tumors, 86 (88.7%) showed immunohistochemical MMR loss and 68 were successfully analyzed by the Promega MSI assay, of which 48 (70.6%) were MSI-high and 20 (29.4%) were MSI-low/microsatellite stable. Seventy-two samples had sufficient DNA for the sequencing-based MSI assay, of which 55 (76.4%) and 61 (84.7%) scored as MSI-high using the 24-marker and 54-marker panels, respectively. The concordance between the MSI assays and immunohistochemistry was 70.6% (p = 0.003), 87.5% (p = 0.039), and 90.3% (p = 1.00) for the Promega assay, the 24-marker assay, and the 54-marker assay, respectively. Of the 11 tumors with retained MMR protein expression, four were MSI-low/MSI-high or MSI-high by the Promega assay or one of the sequencing-based assays. Conclusion: Our results show that Lynch syndrome-associated urothelial cancers frequently had loss of MMR protein expression. The Promega MSI assay was significantly less sensitive, but the 54-marker sequencing-based MSI analysis showed no significant difference compared to immunohistochemistry. Data from this study alongside previous studies, suggest that universal MMR deficiency testing of newly diagnosed urothelial cancers, using immunohistochemistry and/or sequencing-based MSI analysis of sensitive markers, offer a potentially useful approach to identification of Lynch syndrome cases.

Constitutional Microsatellite Instability, Genotype, and Phenotype Correlations in Constitutional Mismatch Repair Deficiency.

BACKGROUND & AIMS: Constitutional mismatch repair deficiency (CMMRD) is a rare recessive childhood cancer predisposition syndrome caused by germline mismatch repair variants. Constitutional microsatellite instability (cMSI) is a CMMRD diagnostic hallmark and may associate with cancer risk. We quantified cMSI in a large CMMRD patient cohort to explore genotype-phenotype correlations using novel MSI markers selected for instability in blood. METHODS: Three CMMRD, 1 Lynch syndrome, and 2 control blood samples were genome sequenced to >120× depth. A pilot cohort of 8 CMMRD and 38 control blood samples and a blinded cohort of 56 CMMRD, 8 suspected CMMRD, 40 Lynch syndrome, and 43 control blood samples were amplicon sequenced to 5000× depth. Sample cMSI score was calculated using a published method comparing microsatellite reference allele frequencies with 80 controls. RESULTS: Thirty-two mononucleotide repeats were selected from blood genome and pilot amplicon sequencing data. cMSI scoring using these MSI markers achieved 100% sensitivity (95% CI, 93.6%-100.0%) and specificity (95% CI 97.9%-100.0%), was reproducible, and was superior to an established tumor MSI marker panel. Lower cMSI scores were found in patients with CMMRD with MSH6 deficiency and patients with at least 1 mismatch repair missense variant, and patients with biallelic truncating/copy number variants had higher scores. cMSI score did not correlate with age at first tumor. CONCLUSIONS: We present an inexpensive and scalable cMSI assay that enhances CMMRD detection relative to existing methods. cMSI score is associated with mismatch repair genotype but not phenotype, suggesting it is not a useful predictor of cancer risk.

Is HLA type a possible cancer risk modifier in Lynch syndrome?

Lynch syndrome (LS) is the most common inherited cancer syndrome. It is inherited via a monoallelic germline variant in one of the DNA mismatch repair (MMR) genes. LS carriers have a broad 30% to 80% risk of developing various malignancies, and more precise, individual risk estimations would be of high clinical value, allowing tailored cancer prevention and surveillance. Due to MMR deficiency, LS cancers are characterized by the accumulation of frameshift mutations leading to highly immunogenic frameshift peptides (FSPs). Thus, immune surveillance is proposed to inhibit the outgrowth of MMR-deficient cell clones. Recent studies have shown that immunoediting during the evolution of MMR-deficient cancers leads to a counter-selection of highly immunogenic antigens. The immunogenicity of FSPs is dependent on the antigen presentation. One crucial factor determining antigen presentation is the HLA genotype. Hence, a LS carrier's HLA genotype plays an important role in the presentation of FSP antigens to the immune system, and may influence the likelihood of progression from precancerous lesions to cancer. To address the challenge of clarifying this possibility including diverse populations with different HLA types, we have established the INDICATE initiative (Individual cancer risk by HLA type, http://indicate-lynch.org/), an international network aiming at a systematic evaluation of the HLA genotype as a possible cancer risk modifier in LS. Here we summarize the current knowledge on the role of HLA type in cancer risk and outline future research directions to delineate possible association in the scenario of LS with genetically defined risk population and highly immunogenic tumors.

Detection of Microsatellite Instability in Colonoscopic Biopsies and Postal Urine Samples from Lynch Syndrome Cancer Patients Using a Multiplex PCR Assay.

Identification of mismatch repair (MMR)-deficient colorectal cancers (CRCs) is recommended for Lynch syndrome (LS) screening, and supports targeting of immune checkpoint inhibitors. Microsatellite instability (MSI) analysis is commonly used to test for MMR deficiency. Testing biopsies prior to tumour resection can inform surgical and therapeutic decisions, but can be limited by DNA quantity. MSI analysis of voided urine could also provide much needed surveillance for genitourinary tract cancers in LS. Here, we reconfigure an existing molecular inversion probe-based MSI and BRAF c.1799T > A assay to a multiplex PCR (mPCR) format, and demonstrate that it can sample >140 unique molecules per marker from <1 ng of DNA and classify CRCs with 96−100% sensitivity and specificity. We also show that it can detect increased MSI within individual and composite CRC biopsies from LS patients, and within preoperative urine cell free DNA (cfDNA) from two LS patients, one with an upper tract urothelial cancer, the other an undiagnosed endometrial cancer. Approximately 60−70% of the urine cfDNAs were tumour-derived. Our results suggest that mPCR sequence-based analysis of MSI and mutation hotspots in CRC biopsies could facilitate presurgery decision making, and could enable postal-based screening for urinary tract and endometrial tumours in LS patients.

How Should We Test for Lynch Syndrome? A Review of Current Guidelines and Future Strategies.

International guidelines for the diagnosis of Lynch syndrome (LS) recommend molecular screening of colorectal cancers (CRCs) to identify patients for germline mismatch repair (MMR) gene testing. As our understanding of the LS phenotype and diagnostic technologies have advanced, there is a need to review these guidelines and new screening opportunities. We discuss the barriers to implementation of current guidelines, as well as guideline limitations, and highlight new technologies and knowledge that may address these. We also discuss alternative screening strategies to increase the rate of LS diagnoses. In particular, the focus of current guidance on CRCs means that approximately half of Lynch-spectrum tumours occurring in unknown male LS carriers, and only one-third in female LS carriers, will trigger testing for LS. There is increasing pressure to expand guidelines to include molecular screening of endometrial cancers, the most frequent cancer in female LS carriers. Furthermore, we collate the evidence to support MMR deficiency testing of other Lynch-spectrum tumours to screen for LS. However, a reliance on tumour tissue limits preoperative testing and, therefore, diagnosis prior to malignancy. The recent successes of functional assays to detect microsatellite instability or MMR deficiency in non-neoplastic tissues suggest that future diagnostic pipelines could become independent of tumour tissue.

Constitutional mismatch repair deficiency is the diagnosis in 0.41% of pathogenic NF1/SPRED1 variant negative children suspected of sporadic neurofibromatosis type 1.

PURPOSE: Biallelic germline mismatch repair (MMR) gene pathogenic variants (PVs) cause constitutional MMR deficiency (CMMRD), a highly penetrant childhood cancer syndrome phenotypically overlapping with neurofibromatosis type 1 (NF1). CMMRD testing in suspected NF1 children without NF1/SPRED1 PVs enables inclusion of CMMRD positives into monitoring programs prior to tumor onset. However, testing is associated with potential harms and the prevalence of CMMRD among these children is unknown. METHODS: Using a simple and scalable microsatellite instability (MSI) assay of non-neoplastic leukocyte DNA to detect CMMRD, we retrospectively screened >700 children suspected of sporadic NF1 but lacking NF1/SPRED1 PVs. RESULTS: For three of seven MSI-positive patients germline MMR gene PVs confirmed the diagnosis of CMMRD. Founder variants NM_000535.5(PMS2):c.736_741delinsTGTGTGTGAAG, prevalent in Europe and North America, and NM_000179.2(MSH6):c.10C>G, affecting 1:400 French Canadians, represented two of five PVs. The prevalence of CMMRD was 3/735 (0.41%, 95% confidence interval [CI]: 0.08-1.19%). CONCLUSION: Our empirical data provide reliable numbers for genetic counseling and confirm previous prevalence estimations, on which Care for CMMRD consortium guidelines are based. These advocate CMMRD testing of preselected patients rather than offering reflex testing to all suspected sporadic NF1 children lacking NF1/SPRED1 PVs. The possibility of founder effects should be considered alongside these testing guidelines.

Sequencing-based microsatellite instability testing using as few as six markers for high-throughput clinical diagnostics.

Microsatellite instability (MSI) testing of colorectal cancers (CRCs) is used to screen for Lynch syndrome (LS), a hereditary cancer-predisposition, and can be used to predict response to immunotherapy. Here, we present a single-molecule molecular inversion probe and sequencing-based MSI assay and demonstrate its clinical validity according to existing guidelines. We amplified 24 microsatellites in multiplex and trained a classifier using 98 CRCs, which accommodates marker specific sensitivities to MSI. Sample classification achieved 100% concordance with the MSI Analysis System v1.2 (Promega) in three independent cohorts, totaling 220 CRCs. Backward-forward stepwise selection was used to identify a 6-marker subset of equal accuracy to the 24-marker panel. Assessment of assay detection limits showed that the 24-marker panel is marginally more robust to sample variables than the 6-marker subset, detecting as little as 3% high levels of MSI DNA in sample mixtures, and requiring a minimum of 10 template molecules to be sequenced per marker for >95% accuracy. BRAF c.1799 mutation analysis was also included to streamline LS testing, with all c.1799T>A variants being correctly identified. The assay, therefore, provides a cheap, robust, automatable, and scalable MSI test with internal quality controls, suitable for clinical cancer diagnostics.

A sensitive and scalable microsatellite instability assay to diagnose constitutional mismatch repair deficiency by sequencing of peripheral blood leukocytes.

Constitutional mismatch repair deficiency (CMMRD) is caused by germline pathogenic variants in both alleles of a mismatch repair gene. Patients have an exceptionally high risk of numerous pediatric malignancies and benefit from surveillance and adjusted treatment. The diversity of its manifestation, and ambiguous genotyping results, particularly from PMS2, can complicate diagnosis and preclude timely patient management. Assessment of low-level microsatellite instability in nonneoplastic tissues can detect CMMRD, but current techniques are laborious or of limited sensitivity. Here, we present a simple, scalable CMMRD diagnostic assay. It uses sequencing and molecular barcodes to detect low-frequency microsatellite variants in peripheral blood leukocytes and classifies samples using variant frequencies. We tested 30 samples from 26 genetically-confirmed CMMRD patients, and samples from 94 controls and 40 Lynch syndrome patients. All samples were correctly classified, except one from a CMMRD patient recovering from aplasia. However, additional samples from this same patient tested positive for CMMRD. The assay also confirmed CMMRD in six suspected patients. The assay is suitable for both rapid CMMRD diagnosis within clinical decision windows and scalable screening of at-risk populations. Its deployment will improve patient care, and better define the prevalence and phenotype of this likely underreported cancer syndrome.

An integrated transcriptional analysis of the developing human retina.

The scarcity of embryonic/foetal material as a resource for direct study means that there is still limited understanding of human retina development. Here, we present an integrated transcriptome analysis combined with immunohistochemistry in human eye and retinal samples from 4 to 19 post-conception weeks. This analysis reveals three developmental windows with specific gene expression patterns that informed the sequential emergence of retinal cell types and enabled identification of stage-specific cellular and biological processes, and transcriptional regulators. Each stage is characterised by a specific set of alternatively spliced transcripts that code for proteins involved in the formation of the photoreceptor connecting cilium, pre-mRNA splicing and epigenetic modifiers. Importantly, our data show that the transition from foetal to adult retina is characterised by a large increase in the percentage of mutually exclusive exons that code for proteins involved in photoreceptor maintenance. The circular RNA population is also defined and shown to increase during retinal development. Collectively, these data increase our understanding of human retinal development and the pre-mRNA splicing process, and help to identify new candidate disease genes.

A novel panel of short mononucleotide repeats linked to informative polymorphisms enabling effective high volume low cost discrimination between mismatch repair deficient and proficient tumours.

Somatic mutations in mononucleotide repeats are commonly used to assess the mismatch repair status of tumours. Current tests focus on repeats with a length above 15bp, which tend to be somatically more unstable than shorter ones. These longer repeats also have a substantially higher PCR error rate, and tests that use capillary electrophoresis for fragment size analysis often require expert interpretation. In this communication, we present a panel of 17 short repeats (length 7-12bp) for sequence-based microsatellite instability (MSI) testing. Using a simple scoring procedure that incorporates the allelic distribution of the mutant repeats, and analysis of two cohort of tumours totalling 209 samples, we show that this panel is able to discriminate between MMR proficient and deficient tumours, even when constitutional DNA is not available. In the training cohort, the method achieved 100% concordance with fragment analysis, while in the testing cohort, 4 discordant samples were observed (corresponding to 97% concordance). Of these, 2 showed discrepancies between fragment analysis and immunohistochemistry and one was reclassified after re-testing using fragment analysis. These results indicate that our approach offers the option of a reliable, scalable routine test for MSI.

Analysis of human ES cell differentiation establishes that the dominant isoforms of the lncRNAs RMST and FIRRE are circular.

BACKGROUND: Circular RNAs (circRNAs) are predominantly derived from protein coding genes, and some can act as microRNA sponges or transcriptional regulators. Changes in circRNA levels have been identified during human development which may be functionally important, but lineage-specific analyses are currently lacking. To address this, we performed RNAseq analysis of human embryonic stem (ES) cells differentiated for 90 days towards 3D laminated retina. RESULTS: A transcriptome-wide increase in circRNA expression, size, and exon count was observed, with circRNA levels reaching a plateau by day 45. Parallel statistical analyses, controlling for sample and locus specific effects, identified 239 circRNAs with expression changes distinct from the transcriptome-wide pattern, but these all also increased in abundance over time. Surprisingly, circRNAs derived from long non-coding RNAs (lncRNAs) were found to account for a significantly larger proportion of transcripts from their loci of origin than circRNAs from coding genes. The most abundant, circRMST:E12-E6, showed a > 100X increase during differentiation accompanied by an isoform switch, and accounts for > 99% of RMST transcripts in many adult tissues. The second most abundant, circFIRRE:E10-E5, accounts for > 98% of FIRRE transcripts in differentiating human ES cells, and is one of 39 FIRRE circRNAs, many of which include multiple unannotated exons. CONCLUSIONS: Our results suggest that during human ES cell differentiation, changes in circRNA levels are primarily globally controlled. They also suggest that RMST and FIRRE, genes with established roles in neurogenesis and topological organisation of chromosomal domains respectively, are processed as circular lncRNAs with only minor linear species.

Interaction between polymorphisms in aspirin metabolic pathways, regular aspirin use and colorectal cancer risk: A case-control study in unselected white European populations.

Regular aspirin use is associated with reduced risk of colorectal cancer (CRC). Variation in aspirin's chemoprevention efficacy has been attributed to the presence of single nucleotide polymorphisms (SNPs). We conducted a meta-analysis using two large population-based case-control datasets, the UK-Leeds Colorectal Cancer Study Group and the NIH-Colon Cancer Family Registry, having a combined total of 3325 cases and 2262 controls. The aim was to assess 42 candidate SNPs in 15 genes whose association with colorectal cancer risk was putatively modified by aspirin use, in the literature. Log odds ratios (ORs) and standard errors were estimated for each dataset separately using logistic regression adjusting for age, sex and study site, and dataset-specific results were combined using random effects meta-analysis. Meta-analysis showed association between SNPs rs6983267, rs11694911 and rs2302615 with CRC risk reduction (All P<0.05). Association for SNP rs6983267 in the CCAT2 gene only was noteworthy after multiple test correction (P = 0.001). Site-specific analysis showed association between SNPs rs1799853 and rs2302615 with reduced colon cancer risk only (P = 0.01 and P = 0.004, respectively), however neither reached significance threshold following multiple test correction. Meta-analysis of SNPs rs2070959 and rs1105879 in UGT1A6 gene showed interaction between aspirin use and CRC risk (Pinteraction = 0.01 and 0.02, respectively); stratification by aspirin use showed an association for decreased CRC risk for aspirin users having a wild-type genotype (rs2070959 OR = 0.77, 95% CI = 0.68-0.86; rs1105879 OR = 0.77 95% CI = 0.69-0.86) compared to variant allele cariers. The direction of the interaction however is in contrast to that published in studies on colorectal adenomas. Both SNPs showed potential site-specific interaction with aspirin use and colon cancer risk only (Pinteraction = 0.006 and 0.008, respectively), with the direction of association similar to that observed for CRC. Additionally, they showed interaction between any non-steroidal anti-inflammatory drugs (including aspirin) use and CRC risk (Pinteraction = 0.01 for both). All gene x environment (GxE) interactions however were not significant after multiple test correction. Candidate gene investigation indicated no evidence of GxE interaction between genetic variants in genes involved in aspirin pathways, regular aspirin use and colorectal cancer risk.

PTESFinder: a computational method to identify post-transcriptional exon shuffling (PTES) events.

BACKGROUND: Transcripts, which have been subject to Post-transcriptional exon shuffling (PTES), have an exon order inconsistent with the underlying genomic sequence. These have been identified in a wide variety of tissues and cell types from many eukaryotes, and are now known to be mostly circular, cytoplasmic, and non-coding. Although there is no uniformly ascribed function, several have been shown to be involved in gene regulation. Accurate identification of these transcripts can, however, be difficult due to artefacts from a wide variety of sources. RESULTS: Here, we present a computational method, PTESFinder, to identify these transcripts from high throughput RNAseq data. Uniquely, it systematically excludes potential artefacts emanating from pseudogenes, segmental duplications, and template switching, and outputs both PTES and canonical exon junction counts to facilitate comparative analyses. In comparison with four existing methods, PTESFinder achieves highest specificity and comparable sensitivity at a variety of read depths. PTESFinder also identifies between 13 % and 41.6 % more structures, compared to publicly available methods recently used to identify human circular RNAs. CONCLUSIONS: With high sensitivity and specificity, user-adjustable filters that target known sources of false positives, and tailored output to facilitate comparison of transcript levels, PTESFinder will facilitate the discovery and analysis of these poorly understood transcripts.

Circular RNA enrichment in platelets is a signature of transcriptome degradation.

In platelets, splicing and translation occur in the absence of a nucleus. However, the integrity and stability of mRNAs derived from megakaryocyte progenitor cells remain poorly quantified on a transcriptome-wide level. As circular RNAs (circRNAs) are resistant to degradation by exonucleases, their abundance relative to linear RNAs can be used as a surrogate marker for mRNA stability in the absence of transcription. Here we show that circRNAs are enriched in human platelets 17- to 188-fold relative to nucleated tissues and 14- to 26-fold relative to samples digested with RNAse R to selectively remove linear RNA. We compare RNAseq read depths inside and outside circRNAs to provide in silico evidence of transcript circularity, show that exons within circRNAs are enriched on average 12.7 times in platelets relative to nucleated tissues and identify 3162 genes significantly enriched for circRNAs, including some where all RNAseq reads appear to be derived from circular molecules. We also confirm that this is a feature of other anucleate cells through transcriptome sequencing of mature erythrocytes, demonstrate that circRNAs are not enriched in cultured megakaryocytes, and demonstrate that linear RNAs decay more rapidly than circRNAs in platelet preparations. Collectively, these results suggest that circulating platelets have lost >90% of their progenitor mRNAs and that translation in platelets occurs against the backdrop of a highly degraded transcriptome. Finally, we find that transcripts previously classified as products of reverse transcriptase template switching are both enriched in platelets and resistant to decay, countering the recent suggestion that up to 50% of rearranged RNAs are artifacts.

Identification of a neuronal transcription factor network involved in medulloblastoma development.

BACKGROUND: Medulloblastomas, the most frequent malignant brain tumours affecting children, comprise at least 4 distinct clinicogenetic subgroups. Aberrant sonic hedgehog (SHH) signalling is observed in approximately 25% of tumours and defines one subgroup. Although alterations in SHH pathway genes (e.g. PTCH1, SUFU) are observed in many of these tumours, high throughput genomic analyses have identified few other recurring mutations. Here, we have mutagenised the Ptch+/- murine tumour model using the Sleeping Beauty transposon system to identify additional genes and pathways involved in SHH subgroup medulloblastoma development. RESULTS: Mutagenesis significantly increased medulloblastoma frequency and identified 17 candidate cancer genes, including orthologs of genes somatically mutated (PTEN, CREBBP) or associated with poor outcome (PTEN, MYT1L) in the human disease. Strikingly, these candidate genes were enriched for transcription factors (p=2x10-5), the majority of which (6/7; Crebbp, Myt1L, Nfia, Nfib, Tead1 and Tgif2) were linked within a single regulatory network enriched for genes associated with a differentiated neuronal phenotype. Furthermore, activity of this network varied significantly between the human subgroups, was associated with metastatic disease, and predicted poor survival specifically within the SHH subgroup of tumours. Igf2, previously implicated in medulloblastoma, was the most differentially expressed gene in murine tumours with network perturbation, and network activity in both mouse and human tumours was characterised by enrichment for multiple gene-sets indicating increased cell proliferation, IGF signalling, MYC target upregulation, and decreased neuronal differentiation. CONCLUSIONS: Collectively, our data support a model of medulloblastoma development in SB-mutagenised Ptch+/- mice which involves disruption of a novel transcription factor network leading to Igf2 upregulation, proliferation of GNPs, and tumour formation. Moreover, our results identify rational therapeutic targets for SHH subgroup tumours, alongside prognostic biomarkers for the identification of poor-risk SHH patients.